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We propose a simple diatomic system trapped inside an optical cavity to control the energy 
flow between two thermal baths. Through the action of the baths the system is driven to a non¬ 
equilibrium steady state. Using the Large Deviation theory we show that the number of photons 
flowing between the two baths is dramatically different depending on the symmetry of the atomic 
states. Here we present a deterministic scheme to prepare symmetric and antisymmetric atomic 
states with the use of external driving fields, thus implementing an atomic control switch for the 
energy flow. 
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INTRODUCTION 

Understanding how energy flows and how to control this flow has recently attracted considerable attention [1]. It 
further prompted the study of transport phenomena in quantum systems including the analysis of fundamental laws 
like Fourier’s Law [2-5] and the second law of thermodynamics [6]. Furthermore, it was shown that quantum effects 
can be used to enhance important thermodynamical processes. Examples include quantum thermal machines [7] and 
solar cells [8-10]. Another exciting field of research is thermal phenomena and design of electronic and thermal devices 
at the nanoscale, where several experiments demonstrated electronic switches at the atomic scale [11-13]. Notably, 
a two-terminal switch to control the charge transport in a junction based on the reversible rearrangement of single 
atoms was shown in Ref. [13]. 

One of the most promising physical platforms to study quantum transport phenomena is coupled optical cavities 
with one or more atoms trapped per site. In order to design different topologies, optical cavities can be linked via 
fibers [14-16] with an effective interaction given by Bose-Hubbard Hamiltonians [17]. These cavities are suitable 
to model transport in systems with next-neighbor interactions with lineal or network topologies [16]. The control 
capability of optical cavities make them a very useful tool to study transport phenomena. Applications include the 
study of quantum transport in different dimensions [18, 19], and the simulation of noise-assisted transport [20]. 

Recently, Buca and Prosen showed that certain symmetries in an open quantum system lead to different nonequi¬ 
librium steady states [21] (see also [22, 23]). These steady states are classified according to the symmetry operator 
spectrum. The multiplicity of steady states leads to different expected currents. This effect can be used to design 
a symmetry-controlled quantum switch to govern the energy current passing through the system by simply selecting 
its initial state [24, 25]. Subsequently, the role of symmetry in energy transfer was studied in both transient [26, 27] 
and steady state scenarios [28]. However, harnessing quantum symmetry to control the energy flux in a steady state 
scenario was considered only in toy models [24]. 

In this paper we propose a feasible design of an atomic symmetry-controlled thermal switch. Our design comprises 
a chain of cavities, coupled, at the two ends, to two different temperature baths, Ti and T 2 , while the middle cavity 
is doped with two laser-driven atoms. The thermal baths drive the system out of equilibrium and it evolves to a 
steady state scenario, where there is a finite energy current flowing from the hot to the cold bath. We show how a 
pair of atoms trapped in the middle cavity can be used to control the energy current through the system up to four 
orders of magnitude. The control is implemented, in a deterministic fashion, by switching between the symmetric and 
antisymmetric atomic state manifolds with the use of a laser driving field. 


RESULTS AND DISCUSSION 

Our system comprises a chain of three optical cavities coupled to thermal baths at temperatures Ti and T 2 , as it is 
shown in Figure (1). The middle cavity is doped with two three-state atoms addressed by external laser fields, while 
the cavity photons hop between neighboring sites at a rate J. The evolution of the system is given by a Markovian 
master equation [29] in the form (h= 1), 
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p = _,^Hp-pH) + Y.^fP+ E (1) 

i=l c={Lr} 

where H = 77'=*''' -p is the rotating-wave approximation (RWA) Hamiltonian accounting for the control atom- 

laser interaction and the hopping of photons, 77 *'°p = J ^a|a 2 -I- a\ar -f H.c.^. Here aij^jr the anhilitation operator 

of photons of the left/middle/right cavity, while Cp = LpL'^ — ^ L, p} are Lindblad superoperators, which describe 
the interaction of the system with the incoherent channels. There are four decay channels associated with the atomic 
states [see Figure (1) right], which are accounted for by Uf* with i = {1,... ,4}, and two decay channels with 
c = {?,r}, associated with the thermal bath coupled to the left and right cavity, respectively. The RWA is taken for 
simplicity, but it does neither affect the multiplicity of steady states nor the control capacity of the system. 

To realize a controlled deterministic switching between atomic states with different symmetry we adopt a scheme, 
where the quantum state preparation is heralded by macroscopic quantum jumps [30, 31]. We assume that two 
identical three-state A atoms are coupled to the middle cavity’s mode 02 on the |0) O |e) transition with a coupling 
strength g. Furthermore, the |1) O |e) and jO) <->■ jl) transitions are driven by two external laser fields with Rabi 
frequencies Uq and Ui, respectively. If the excited states |e) are far off-resonant, i.e. Hi < g, F, Ug where A and 
r are the detuning and decay rate of the excited states, then they can be adiabatically eliminated. Thus, the state 
space of the diatomic system is reduced to contain only four states, three of which jOO), |11), and |s) = (|01)-|-|10))/-\/2, 
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FIG. 1. Left: A schematic of a symmetry-controlled switch to control the current between two thermal baths, Ti and T 2 . Right: 
Energy diagram of the symmetric and antisymmetric atomic subspaces, which modify the energy current between the baths up 
to four orders of magnitude, and the possible transitions between them. 


are symmetric, while |a) = (|01) — |10))/-\/2 is antisymmetric, as displayed in Figure (1) right. The system is then 
equivalent to a four level system. Then, the control Hamiltonian is given by (see Methods section) 

(|00)(s| + |s)(ll| + H.c.) + g' (|00)(s|4 + |s)(ll|4 + H.c) + A' (|00)(00| - |11)(11|), (2) 

with A' = -^ 4^2 - and 5 ' = 

The transitions between the symmetric and antisymmetric atomic states occur as a result of two decay channels, 

Lf = VK [|00)(a| - |a)(ll|], Tf = ^ [b)(a| + |a)(s|], (3) 

while there are two more decay channels which link the symmetric states, = •4]4[|00)(s| -I- |s)(ll|] and = 

4^ [|Q)(a| + |s)('S| + 2|11)(11|], see Figure (1) right. Here, the decay rates are modihed by the external driving helds 

r 

and are explicitly given by = " 4^2 ° , with Fq -I- Fi = F, [30, 31]. Note that, there is mixing between the 

symmetric and antisymmetric atomic subspaces only when Hq 7 ^ 0 . 

The presence of the baths leads to the creation and destruction of cavity photons in the left and right cavities, 
which we model by the following incoherent channels, 

= \/7th(ni -b 1)4: 

T ‘2 = \/7th(n2 + 1) 4: (4) 

where subscripts 1 and 2 correspond to creation and destruction, respectively, 7 th is the interaction rate, ni( 2 ) = 
l/[exp(tLi/(fcBT 4 ( 2 ))) — 1 ] is the temperature-dependent mean excitation photon number at the resonance frequency 
in the respective bosonic bath, and fee is the Boltzmann’s constant [29]. 

To analise the transport properties of our system we use a Large Deviation approach [24, 32-34] . The key 
elements of this technique are the large deviation functions 0{s) and G{q) that account for the statistic of photons 
interchanged between the system and one of the baths. These functions determine the values for long times of both 
the moment-generating function, Zgit) ^ exp[t 0 (s], and the probability of having a certain value q of the current. 
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FIG. 2. Large deviation functions 6{s) (main) and G{q) (inset). In both plots the blue line represents the state with the laser 
on (Hi = O.OOSg, Ho = (?) and the red line with the laser off (Hi = Ho = 0). The remaining parameters are A = 75g, To = 
Fi, r = g, J = ni = 0.005, n 2 = 10“®. The dashed line in G{q) sketches the non-convex regime that cannot be 

calculated from 9{s). 


Pq{t) ~ exp [t0(s)] (see Methods). The moments of the current distribution can be calculated directly from the large 
deviation function 9{s) by the expression (g“) = — ^ 

S — {) 

In Figure 2 the moment-generator LDF 0{s) is displayed for the cases when the lasers are driving the atoms 
(flo.i 7 ^ 0) and when they are off (flo.i = 0). When the lasers are off there are two different steady states in the 
system. This is reflected in the kink of 0{s) at s = 0. The non-analyticity of 0{s) can be interpreted as a dynamical 
phase transition, and it is a consequence of the existence of more than one steady state with different activity. If 
the atoms are in one of the symmetric states {| 00 ), | 11 ), |s)}, then they facilitate the transfer of photons to the cold 
reservoir Ti. This is described by 6{s) in the region with s < 0. However, if the atoms are in the antisymmetric state 
|a), which is a dark state for the system, then they cannot interact with the cavity photons, thereby reducing the 
energy transfer through the cavities. This scenario is depicted in the LDF 0{s) for the regime s > 0. 

The multiplicity of the steady states is a consequence of the symmetries of the system [24]. If Oq = 0 (laser off) the 
unitary operator tt that interchanges the states of the two atoms (TrjOO) = jOO), 7 r|s) = |s), Trjll) = jll), 7 r|a) = —|a)) 
commutes with both the Hamiltonian and all of the Lindblad operators, [tt, 7J] = [tt, Lf*] = [ 7 r,L[:*'] = 0, V(i, c). 
This leads to a degeneracy of the Liouvillian operator and to multiple steady states, which can be labeled by the 
eigenvalues of the operator tt [21]. When the laser is on, meaning that Dg 7 ^ 0, the degeneracy is broken, because the 
operator Lf- from Eq. (3) does not commute with tt because [tt, Lf*] = 2Lf‘ is equal to zero if and only if Lf* = 0, 
that corresponds with Hg = 0. In this case there is only one steady state that mixes the symmetric and antisymmetric 
manifolds of the atomic states. Because of this, there is no dynamical phase transition and 9{s) and its derivative are 
analytical for all values of s, including zero, as it is displayed in Figure 2. 

We recover the large deviation function G{q) from 9{s) by a Legendre transformation. The non-analyticity of the 
derivative of 9{s) leads to a non-convex regime in G{q) that cannot be directly inferred by the Legendre transformation, 
as it is sketched in the inset of Figure 2. The probability of having Q events in a long time t can be calculated as 
PQ{t) = exp[tG(Q/t)]. 

To compare the efficiency of the symmetric and antisymmetric steady states we calculate the mean number of 

. If there are more 


photons transferred to the reservoir per unit time (g) directly from 9{s) according to (g) = — 


s=0 


than one steady states, the mean flux for the more and less active steady states are given by (g)max = — 


d0{s) 

ds 


(g)min = - 


dS{s) 


ds 


s—>-0 
s<0 


and 


s—>-0 
s>0 


, respectively. We calculate the ratio between the maximum and minimum flux a = (g)max/(g)r 


In Figure 3 we plot a as a function of the hopping parameter J and we find that a is strongly dependent on the value 
of J, up to three orders of magnitude. The control capacity of the system increases dramatically when J decreases. 
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FIG. 3. Left, blue: Ratio between maximum and minimum flux as a function of the hopping parameter J. We have assumed no 
driving of the atoms (fti = Hq = 0). Left, red: Ratio between the mean length of a dark period and the mean time between the 
absorption of photons in a bright period (Hi = 0.0025(7, flo = 5 )- Right: Total flux for the bright state (fli = 0.00253, flo = <?) 
as a function of the hopping parameter J. The flux is expressed in atomic units and it corresponds to ~ 10^^ photons per 
second. In all plots the remaining parameters are Lq = Ti, T = 3 , ni = 0.1, n 2 = 0. 


and vice versa. The reason is that J modifies the coupling of the atoms to the cavity modes - the smaller J is the 
stronger the atoms are coupled to the photons- which leads to an increase in the control efficiency of the current. 

To control the switch we need to be able to distinguish between the symmetric and antisymmetric (dark) states. 
For our parameter selection the dark state current is close to zero, and the probability of photon absorption for this 
state is negligible. This means that if one photon is absorbed by the cold reservoir we are in a bright state with a 
very high probability. To be able to detect it we need to ensure that the time between absorption of photons in the 
bright state is smaller than the mean length of the dark period. The probability of transition from a dark into a light 
period is Tq + ^T'^, and the mean length of a dark period equals [31], 


Td 


1 8A2 

2To+Tili[- 


(5) 


The mean time between photon absorption can be calculated as the inverse of the mean flux Tc = l/{q). As the mean 
flux is an average of the bright and dark fluxes we can consider it as a lower bound of the flux when the atoms are in 
the bright manifold. In Figure 3 (left) we can see that the ratio Tc/Tu increases with the hopping between cavities J. 
There is a competition between the capacity to distinguish the bright and dark periods that require a high current, 
and the ratio between the bright and dark flux that is higher for smaller fluxes. Nonetheless, for a broad range of the 
parameters the time between photon absorption in the bright state is long enough, compared to the mean length of 
the dark period and the photon flux of the bright state is more than two orders of magnitude larger than the photon 
flux in the dark state. 

The maximum flux, (g)maxj is displayed in Figure 3 (right) as a function of the cavities hopping parameter J. Due 
to the small value of the coupling parameter J, the calculated values for (Q')max are also small but, notably, they are 
measurable. To assess this, we performed a simple calculation and obtained a ((7)max 10“® (in atomic units) that 

corresponds to a flux of the order of 10^^ photons per seconds. Given such light intensity for the bright period, it 
will be easily distinguished from the dark period when the flux is several orders of magnitude smaller. The proposed 
design has several potential applications. It can be used as an energy current switch, an atomic memory and also 
an entangling scheme due to the indirect measurement of the atomic state through the measurement of the current. 
This is particularly useful for memory and entanglement generation, as the antisymmetric state of the atoms is a 
maximally entangled state. Thus, by monitoring the energy current we can control the separability properties of the 
atoms. The readout of the atomic state can be performed without directly interacting with the atoms themselves, as 
their state affects the macroscopic light flux. The proposed system can also be extended to realize better control by 
placing more pairs of atoms in the cavities. By connecting these cavities with optical fibers we can create networks 
of cavities [16], each of them acting as a switch, thereby increasing the flux control. 

In this paper we have presented a realistic model of an atomic thermal switch based on coupled optical cavities. We 
showed that by placing two laser-driven atoms in an optical cavity it is possible to control the energy current trough 
the system when it is driven out of equilibrium by the action of two thermal baths at different temperatures. 
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The control capacity of the system depends on the system parameters, increasing when the hopping between cavities 
decreases. For small values of the hopping parameter between the cavities the current of the bright state can be three 
orders of magnitude bigger than the current of the dark state. This can be used also as a quantum memory as the 
state of the atoms that can be measured without distorting the atoms, just by measuring the current between the 
full system and the thermal bath. The importance of this design relies on the simplicity of manipulating an internal 
degree of freedom of a system in comparison with modifying the system parameters. Similar symmetry-controlled 
switches can be designed for many different quantum setups such as quantum dots, optical lattices and trapped ions. 


METHODS 

Calculation of the Master Equation.- The symmetry control of our system is given by a pair of A atoms at the 
second cavity driven by a laser. This scheme was first proposed to create entangled states by macroscopic quantum 
jumps [30, 31]. The Hamiltonian of the system can be decomposed in the way H = ^ctri _|_ ^hop^ were is the 
part of the Hamiltonian that models the interaction between the atoms and the laser field and H^°p is the photon 
hopping between cavities. The atoms are trapped in the cavity in a way such that the effect of the laser is the same 
in both of them. The coupling strength of the 0 — 2 transition to the cavity field is g, and the laser Rabi frequencies 
for the jo) — | 1 ) and | 1 ) — | 2 ) transitions are Hq and Hi, 


2 2 

= E (^ + H.c) + A|2),(2|, ( 6 ) 

where a\ is the bosonic ladder operator that creates and destroys a photon at cavity 2. It is assumed that Hq < 
g, r. Hi ^ A. As the detunning A is bigger than the other parameters of the system the excited atomic states can be 
adiabatically eliminated. The system is then equivalent to a four level system. As both atoms are equally affected by 
the laser, it is convenient to use the Bell basis for describing the state of the two atoms, with |s) = (jOl) -I- |10))/-\/2 
and |a) = (jOl) — |10))/-\/2. The control Hamiltonian becomes 


^ ^ ^ ^ g, (|00)(s|4 + |s)(ll|4 + H.c) + A' (|00)(00| - |11)(11|), (7) 


with A' = —^0202 — 1^, and g' 
are 


vAa ■ 


The Lindblad operators that correspond to the atoms incoherent channels 


Lf = v/T^ [|00)(a| - |a)(ll|], Lf = [|s)(a| + |a)(s|] , 

Lf = [|00)(s| + |s)(ll|], Lf = ^ [|a)(a| + |s)(s| + 2|11)(11|] (8) 

where T' 

Large deviation approach.- To analyse the statistics of the current flowing between the system and the thermal 
reservoirs we use a full-counting statistics method [24, 32-34]. For simplicity, we focus only on the current between 
the system and the T2 reservoir. We first introduce the reduced density matrix pgit), which is the projection of the 
full density matrix onto the subspace with Q photons interchanged between the system and the T 2 bath. It is not 
normalised and the probability of having a certain value of the current Q in a period of time t is given by Pq (t) = 
Tr[pQ(t)]. This probability scales for long times following a large deviation principle, Pgit) ~ exp[tG'(q)] {g = Q/t), 
where G{q) is the LDF of the current. 

The direct calculation of G{q) is complicated, but it becomes simpler after a change of ensemble. We make 
a Laplace transform on pq(<) and obtain a new reduced density matrix ps = pQ(t) exp(—sQ), where we have 
introduced a counting field s, that is the conjugate field of the number of photons Q. The fields Q and s are dynamical 
variables with the same relation as the thermodynamical variables pressure and volume. This procedure is formally 
equivalent to the thermodynamic transformation between the canonical and microcanonical potentials. The utility 
of this transformation comes from the fact that the evolution of the reduced density matrix pg unravels into a set of 
equations in the form 

4 

Ps = —i {Hps — PsH) + ^ Gf"Pg -\- Cffpg + G^j^Pg + -I- G^^ps = WgPg, 

i=l 


( 9 ) 
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where = e~^Ll^pL^^ — ^ p^ and C^^p = e^Ll^pL^^ — ^ {L^^LI^, p}, and they account for the introduction 

of the counting field s in the right current. 

The Laplace transform fulfills a large deviation principle in the form Zs{t) ~ exp[t0(s)], with 9{s) being a large 
deviation function related to G{q) by a Legendre-type transform 0{s) = maXg[G'(g) — sq\. Here we have introduced 
Zs{t), which is the generating function of the current momenta. The LDF 9{s) corresponds to the eigenvalue of the 
superoperator Wg with the largest real part. The existence of more than one steady states with different currents 
leads to a non-analytic behavior of the LDF 9{s) at s = 0 [24]. 

To calculate the statistical properties of the current we need to diagonalize the superoperator Wg. However, this 
is a numerically intractable problem as the dimension of our system is infinite due to the infinite number of photon 
modes. To overcome this challenge, we adopt a low photon number approximation, where the number of photons 
in the system is small, which can be realized if the temperatures of both baths are very small and the frequency 
of the laser Dq is also small. In this regime we assume that there can be at most one photon in the total system. 
Under this assumption, the Hilbert space required to describe the system is equivalent to the Hilbert space of a four 
qubit-system. Note that, we make the low photon number approximation only for numerical tractability. Our scheme 
will operate as a current switch independent of the number of photons in the system. 
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